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Abstract 

A selfconsistent system of equations describing evolution of linear 
spherically symmetric perturbations in Fridmann world by an arbitrary 
equation of state is obtained. The perturbations singular part correspond- 
ing to massive particle-like source is extracted, an evolutionary equation 
for mass of this source is obtained and exactly solved. An exact solution of 
evolutionary equations of state for perturbations by an arbitrary equation 
of state is constructed. 



1 Introduction 

In a number of papers of one of the Authors collaboratively with A. A. Popov 
Ref.[l]-[3] a theory of spherical perturbations of the Friedmann world in con- 
nection with the necessity of developing the relativistic kinetic theory consid- 
ering gravitational interactions was constructed. The procedure applied for 
the kinetic equations derivation of averaging local fluctuations of the gravita- 
tional field pointed out by one of the Authors an interesting fact: the average 
quadrature fluctuations of the gravitational field play a role of tensor of energy- 
momentum of an ideal fiuid with extremely- hard equation of state [3]. In the 
paper of one of the Authors collaboratively with A. A. Popov exact solutions of 
the equations for the spherically symmetric perturbations of the ultrarelativistic 
Friedmann world with an arbitrary curve index [4] were obtained. However, the 
analysis of the obtained solutions and the calculation of the average quadrature 
corrections to the Friedmann metrics were not fulfilled. In connection with the 
problem of dark energy and dark materia in cosmology the task of possible al- 
teration of macroscopic equation of state of the Friedmann universe by the local 
gravitational interactions becomes actual, as far as the potential energy on the 
gravitational interactions is negative, it can correspond to a negative macro- 
scopic pressure. In the present work we start investigating this possibility by 
fulfilling subsequent maintaining a dynamic theory of the Friedmann universe 
subject to the local gravitational perturbations. 
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2 Spherically symmetric space-time 

2.1 Space symmetry and energy- momentum tensor 

Let us study the space-time with spherical symmetry which metrics in the 
isotropic frame of referenc^ll (r,^, 1^,77), where 77 is the time coordinate and 
r is a derivative one, can be written in the form 

ds^ = e'^drj^ ~ e^[dr^ + r^{d0^ + sin^ 0d(p^)], (1) 

where 

A = A(r,77); v = v{r,ri) 

are arbitrary scalar functions of their arguments. Isotropic coordinates are 
convenient because within these coordinates the three-dimension space metrics 
takes an explicitly conform plane form. As it is known the metrics ([T]) admits a 
group of G3 rotations with the Killing's vectors (see, for example Ref. [7]) 

= (0, siniy9, 6'cos(/5, 0), = (0, — cosip.O sinip, 0), 

r = (0,0,1,0), ' 

3 

Therefore in consequence of the Einstein equations the space-time symmetry is 
taken after by the energy-impulse tensor 

Lr'==0, (a^M), (3) 

where L is the Lie's derivative in the direction of ^ (see, for example Ref. [7]) 
LA' = A^i!" + A^S^j^ (4) 

the other tensor derivatives are processed in the same way as the Lie's derivative. 
Thus under the conditions of spherical symmetry of space-time the energy- 
impulse tensor, as it is known, takes the form of the energy-impulse tensor of 
an ideal isotropic fluid 

T'^ = {e+p)u'u^ -pg'^, (5) 

Where the scalars £(r, 77) and p{r, 77) are an energy density and a fluid pressure 
respectively and zi* is a singular time-like vector of this fluid dynamic velocity 

giku'u'' = 1, (6) 

whereas 

= {u-{r,'n),(),0,u\r,7^)). (7) 



^That is in the frame of reference in which the metrics of three-dimention space takes a 
conform plane form 
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Assuming 

v!' = vvPe 2 , < 1, 
where v(r^r\) is a radial three-dimension velocity of fluid, we get from (|6]) and 



m'' = e 2 



e + V p 
1-^2 ' 



(8) 



1_^2 ' ^2 - ^3 - P- 



Thus the following algebraic relations are valid 



(9) 



T4-ri = (£ + p)i 



(10) 



2.2 Einstein Equations 

Nontrivial Einstein equations relatively the metrics ^ have the form (see, for 
example Ref. [6]jl: 



ev + p 



(11) 



2(A" + //") + J^''+-(A' + z.') 
r 



Snp (= -SttTI), (12) 



-e - ( A" + i A'^ +^A' ) + ^e- A^ = 8.i±^ (= 8nTt); 
4 r / 4 1 - 1;2 

ie-\2A' - z.'A) - 8^(£ + p)e^''-^^/^-^ (- SttTI), 
2 1 — 



(13) 



(14) 



^In order to obtain these equations in the formulae of Ref. [6] it is necessary to assume 
/i = A + 21nr . 
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where /' is a derivative from the function / by the radial variable r and / is 
a derivative by the time variable 77 and the universal set of units is assumed 
everywhere in which G ~ c ^ h ~ 1. Subtracting from the both parts of Eq. 
PT|) the corresponding parts of Eq. (IT^ subject to ^ we get the consequence 



\ +AV' - i v''' +-{>! + v') - (A" + v") 
I 2 r 



^S^le+p)^—^. (15) 

3 Background space-time 

Whithin the isotropic spherical coordinates the nonperturbed gravitational field 
is described by the metrics of homogeneous isotropic universe 

ds^ = a^M ( dif - -^[dr^ + r^(de^ + sin^ ed^^)]] , (16) 
V P W J 

where 

p{r)^l + hr\ (17) 

And the curvature index k — is for the space-plane universe, k = ±1 is 
for the universe with positive and negative curvature of thee-dimensional space 
respectively. Whereas r and are dimensionless variables and the scale factor 
0(77) has length dimension. 

Thus in a nonperturbed state the introduced scalar metrics functions A and 
u are equal to 

lyo^lna^iv); Aq = In (^^^ , (18) 

Whereas in consequence of time space homogeneity 

P^Poiv)'-, e==£o('7); vo^O. (19) 

Substituting Eqs. ([TSl) and (fT9|) into the Einstein Eqs. (fTT |) -(fT5 l) we obtain a 
set of equations describing the Friedmann universe dynamics 



1 fd^ ,\ 8tt 

2 



3 ^; (20) 



1 /^a d^ 



,2 



2 T + A: = -87rp. (21) 



As it is known the second of these equations ((2T|) can be substituted by the 
algebraic-differential consequence of 1 and 2 (see, for example Ref. [6]) 

£ + 3-(e+p) = 0. (22) 
a 
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Here it is convenient to pass on from the time variable 77 to the physical time t 
by the formula 

a{r])dri = dt, f a[ri)dr]\ (23) 

whereas 

1=4- / = «/.. (24) 

Then the independent Einstein equations can be written in a more compact 
form (see, for example Ref. [6]) 



1 /■2 ,N StT 



2(d?+/c) = ^£; (25) 



et+3^(e+p) = 0, (26) 

If we know the equation of state, that is the relation of the form 

P = p{e), (27) 
then Eq. is integrated in quadratures 

de 

-3 In a = f + Const. (28) 

Substituting the solution into Eq. (I^Hl) we obtain a closed differential 
equation of order 1 relatively £(77). In the case of barotropic equation of stat^ 

p = (29) 
Eq. (IM)) is easily integrated 

e = cia-3(«+i), (30) 
and Eq. ([^0]) is integrated in quadratures 

= C2?7, (31) 



^a2-3(«+i) 



where ci and C2 are arbitrary constants. The pointed out equations are inte- 
grated in the elementary functions for the early Universe (t — 0). As it is known 
in this case the behavior of the solutions does not depend on the curvature k 
(see, for example Ref. [6]) and does not differ from the behavior of the solutions 
for a space plane Universe {k = 0) 

a = air]^ ; e = cia^^'-^^''\~^^ , k + 1 ^ 0, (32) 



^We call the reader's attention to the difference , A; is a curvature index and k (kappa) is 
a bartrop coefficient . 
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Where the constants ai and ci are connected by the relation 



3k + 1 



ai = |^(3K + l)^^j , l + 3«^0. (33) 
In doing so let us reduce a relation between the variables t and rj from Eq. (|23p 

3k + 1 3(1 + k) 



Subject to (1321) and dSJ) let us reduce from 1^ 
1 



(35) 



27r(K+l)2t2 

Note the solution of the Einstein equations by k = —1/3 is specific 

/M — 1 

a = C2eV 3 cii)^ (-^ ^ _ ^ g (-cx),+oo). 

o 

However this property is a coordinate only. Really, by passing on from the time 
variable 77 to the physical time t 



we obtain a ^ t and the energy dense formula psp in which it is necessary to 
substitute k = — 1/3 only. By k = — 1 we get from ([5^ the so-called inflational 
solution 

a = --^, f=-ln77; 

At 3A2 /r»y^\ 

a = aie ; e = — — = const. (36) 
on 



The solutions (|33| corresponding to the values k < — 1 describe the so-called 
dark materia. 



4 Linear spherically-symmetric perturbations of 
Friedmann space- time 

4.1 Equations for spherically-symmetric perturbations 

Let us consider now the small spherically-symmetric perturbations of isotropic 
cosmological solution assuming 

(37) 

P = Po{v)+-£ Se; e ^ eo{f]) + Se, 
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Where the scalar functions SX{r, 77), Sv(r, 77), Se{r, 77)and 77) will be assumed 
small of order 1 by smallness. Substituting (l37l) into Eq. (fT5|) in the first 
approximation by smallness of the perturbations SX, Sv, de and v we obtain one 
closed equation relatively the function SX + Sv 



d_ ( p{r) 
dr 



= 0. (38) 



Integrating Eq. (pS)) we get (see also Ref.[4]) 

r Ci(ry)+C2(?7)r2, fc = 0; 
X + u=\ (39) 

where Ci(77)and 6*2(7?) are arbitrary functions. 

Further we shall search only the solutions of the Einstein equations of class, 
which out of some sphere coincide with homogeneous isotropic nonperturbed 
solution 

Such solutions correspond to the retarded solutions of hyperbolic type equations. 
The physical sense of the solutions we shall discuss later. Then according to 
([39)) it should take place (see Ref. [4]) 

5X + 5u^Q,^ 6X^~5v. (41) 

In the paper [4] spherically symmetric perturbations are studied in the ultrare- 
lativistic universe (k = 1/3) only, however in doing that solutions of linearized 
Einstein equations for all types of Friedmann universe were obtained. In the 
present paper we will restrict ourselves with the case of the space-plane universe 
(fc = 0), but the barotrop coefficient k will be considered to be arbitrary. Thus, 
subject to the background Einstein equations ([20| -([2T |) . which in the case of 
fc = have the consequence 

2- = 4(1 -3^), (42) 

we get a closed set of three differential Einstein equations linearized around 
the background solution (|16|) . relatively the three unknown variables Si'{r,ri), 
fc(r, ?7)andw(7-, 77) 

• 2 

Si) + 3Siy- - 3KSiy% = Sna^Sp; (43) 
a 

36i>- +3Siy% ~ ^-^r^-^5v = -STra^fc; (44) 
a or or 
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1 d 

a(5i^' = -87reo(l + K)f ■ (45) 

a-^ or] 

The latter of this set of equations (|45l) is a definition of the radial velocity v{r, rj) . 
One of the equations (|43|) and (j44|) determines the density energy perturbation 
Ss(r, rf). 

4.2 Derivation of particlelike solutions 

Further assuming the investigation of particle-like solutions of perturbation 
equations also, let us study canonic equations of motion of gravitating clas- 
sical point particle in the gravitational field to which the i5-like energy density 
corresponds. As a result of the two competitive processes - the accretion of ma- 
terial environment and the reverse process - evaporation of substance the mass 
of the classic point particle in the material environment cannot be constant. 
Therefore let us write the Hamilton invariant function of massive particle in the 
fornfl 



H(x, P) = y/g^'^P^Pk - m, (=0), 

where m = m{s) is a scalar function. From (j45p we get the normalization 
ratio 

{P,P)^m^{s). (46) 
The relativistic canonical equations of particle motion take the form 
dx' _ dH dP, _ dH 

'ds""dPi' lb'~'dx^' ^ ' 

From the first couple of the canonical equations subject to the normalization 
ratio we obtain 

dx"^ du" du^ , 

-^ = — ^9^k-, — i- = l- (48) 
as ni ds as 

The second couple of the canonical equations of motion gives the Lagrangrian 
equations of classical massive particle of variable mass 

In spherically symmetric metrics the motion equations solution (pS)) which does 
not break the spherical symmetry, is a time line which corresponds to the particle 
rest state at the origin of the coordinates 

r = 0, x'' = 7?, (50) 



''See the details in Ref. [8]. 
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in doing so the mass at rest can be an arbitrary function of the coordinate time 

m = m{ri). (51) 

Let us write down the density energy corresponding to its singular part in in- 
variant form 

fe„, = m{ri)5{r), (52) 

where 5{r) is Dirac invariant (5-function in spherical coordinates perceived in 
the sense of integral ratio 

! d^VS^ix) ^a^Jdn / r^dr6{r) = Aira^ / S{rydr = 1, (53) 
Jo Jo 

so that 

Jd'^VSe^.^mirj). (54) 

Temporally abandoning the time derivatives by rj in the left part of Eq. (j43p 
we get then the following equation for the singular part corresponding to the 
singular part of density 

r —dv — STTa mlii^ovr). (55) 
or ar 

Multiplying the both parts ([55]) by ar'^dr and integrating, then integrating by 
parts in the equation left part and at the same time assuming 

limr^^^O (56) 

r^o dr 

we obtain 

ar'^—^ = 2m,{r]). (57) 
or 

Integrating this equation one more time we get 

2m(ri) 
Ov = . 

ar 

Thus the ratio similar to the known one subject to redetermination of the in- 
variant (5-function (j53p takes place 

;i|:^'|:(-7)=4™'"^^W. (58) 

Therefore in order to extract the particle- like singular part of the solution further 
on it is convenient to introduce a new field function tp{r,ri) such as [9] 

5^ = .5X = ^thjll^^ ^ (59) (59) 
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At that according to ([5^ the relation must be fulfihed 



lim — 

r— >o r 



< oo. (60) 



Extracting the singular part of energy density in the right part of Eq. (|43p , sub- 
stituting the function Si/ in the form (|59|) into Eqs. (j42 |) - (|44l) and excluding the 
singular part subject to the relations ([551) and ([50]) we get the linear equations 
set relatively the function $ and perturbations of energy density and velocity 

$ + - -(1 + = 47rra^Kfc, (61) 

a 2 a- 

3^4 - V" = -4Wfc, (62) 

= -47rra^ 1 + k koi'- (63) 

ar r 



Multiplying Eq. ([62]) by k and adding its both parts to the corresponding parts 
of Eq. ([61]) we obtain a closed equation 



3 

$ + -(1 + 3k)$ - -(1 + k)^<I> - K^j" = 0. (64) 
a 2 



Then assuming according to ([59]) 

^(r, 77) = '0(^1 ~ »77,(77) (65) 



And dividing the variables in Eq. ([64p . we get two equations for the functions 
m{ri) and ipir, rf) 

a 3 

m H — (1 + 3K)r7i — -(1 + K)^m — QM; (66) 
a 2 

o • 2 

V; + ^(1 + 3k)?/> - -(1 + - K^j" = -9(77), (67) 

where 6(77) is an arbitrary function of its argument. 



4.3 Basic theorem 

Further let m = M(9, 77) be the private Eq. ([5S[) solution corresponding to the 
given function 6(77). Then in consequence of Eq. (|67p linearity the function 
ipi — — M(9,77).is a private solution of this equation. Then Eq. ([67| general 
solution can be written in the form 

^-(7^,7/) =7/.o(r,77)-M(e,77), (68) 



10 



where vE'(r, 77) is general solution of the corresponding homogeneous equation 



d 3 d^ 

\E' + -(1 + 3k)*--(1 + k)^1'-k*" = 0. (69) 
a 2 

Further on in consequence of the Eq. (j66p linearity its general solution 

is sum of the corresponding homogeneous equation ^{r/) and private solution of 

inhomogeneous one 

m{7]):^ H{r])+M{e,7]). (70) 

But then 

$(r, 77) ^ ip{r, 77) - m(T?) = *(r, 77) - ^(77), (71) 

where /i(r/) is the homogeneous equation general solution 

d 3 d^ 

/i + -(l + 3K)/i--(l + K)^/^ = (72) 

The rest equations of the set (pT|) - (p5)) describe the evolution of nonsingular 
part of energy density and perturbations velocity 

9 



97- 



= -ATrra'\l + k)£ov. (74) 



Thus we proved the theorem 

Theorem. Linear spherically symmetric perturbations of Friedmann metrics 
are described by a set of two independent linear homogeneous equations (j69p 
and (j72[) relatively two functions ii{ri) and 'ijj{r, rj) which are nonsingular at the 
origin of the coordinates. Spherically symmetric perturbations of energy density 
and velocity are determined through metrics perturbations by the relations (1731) - 
(d . 

At K > the homogeneous Eq. ()69|) is hyperbolic, at k, < it is elliptic, at k — 
this equation coincides with Eq. (|72[) . 



5 Evolutional equations for perturbations at con- 
stant barotrop coefficient 

5.1 Mass evolution of particle-like source 

Let us study the cosmological mass evolution of particle-like source. Passing on 
from the variable rj to the variable 0(77) subject to the relation (|42|) in the mass 
evolution equation ()72|) let us reduce it to the form 



d^/i ^ 3 1 + K d/i ^(1-)- ) 
da^ 2 a da 2 a' 
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The general solution of this equation can be easily obtained 



^l^C^a + C-a-^^^+''\ (76) 

where C+ and C_ are arbitrary constants. The C+ coefficient term in this 
solution corresponds to accretion processes, the C_ coefficient one corresponds 
to evaporation processes. Substituting in Eg . ([751) the scale factor relation to 
the time variable 77 ([5^ and using the relation (IM)) between the time variable rj 
and the physical time t we obtain in an explicit form the law of mass evolution 
of particle-like source 



/i = C+i3(i+'') +C-t~\ (i + k)7^0, (77) 

Where C+ and C_ are new arbitrary constants. From the solutions ([75| and 
(|77|) it is seen C_ = C_ corresponds to the final particle mass at the moment 
of time i = 0. The solutions (|76|) - ([77| summarize the solutions obtained in the 
previous papers [3]- [4] for the two private values of adiabatic curve coefficient 
K = and K — l/3which correspond to nonrelativistic and ultrarelativistic 
equation of state correspondingly. In the pointed out cases we get from ([TT)) 

^jL^ C+tl +C-t-^, K = 0, (78) 

li = C+t\ +C-t-\ n=\, (79) 

Let us study a numerical example. Let the particle-like source mass equal to 
Planck mass at the moment of time t — tpi then at the contemporary moment 
of time t ^ lif'^tpi according to ([77)1 the contemporary mass of "particle" 
varies within the limits 10~^^Mq at k = 1 till 10^ Af© at k = 0. At negative 
values of the adiabate curve index the perturbation mass rapidly increases and 
at K « —0, 5 in order of value it is compared with the visible universe mass (see 
Fig.l). 

In the case of the inflation solution the evolutional equation (|72p takes the 
form 

2 

il = — /i = 0. (80) 
r; 

Solving ([50)) we obtain 

m = C+ + — . (81) 

■n 

Now using the relation of the time constant 77 to the physical time t for the case 
of inflation ([55)1 for this case finally we get 

/i = C+ C-e-^K (82) 



12 



In this case C_ = 



corresponds to the final mass at i — > — cxiagain. 



40- 



20 



10 




Fig.l. Relationship of contemporary mass value of particle- like source in Fried- 
mann world which had Planck mass at Planck moment of time to barotrop 
coefhcient k. On the ordinate axis the mass logarithm values are marked off in 
Solar mass units Mq k, 2 • 10'^^ g. 



5.2 Evolutional Equation for nonsingular mode of pertur- 
bations 



From the general nonperturbed solution p2[) we obtain useful ratios 
d 2 1 h 2(1 -3k) 1 



3k + 1 77 ' 



(3k + 1)2 ,y2 



(1 + K^O), 



(83) 



using which we reduce Eq. (IS^ for nonsingular mode of perturbations to the 
form 



■■ 2. 6(1 + k) * ^ 
?/ (1 + Any r]^ 



(84) 



In the case of (1 + k) = Eq. (|74|) is no longer an equation for defining a 
radial velocity of perturbations, but becomes a differential equation relatively 
the function ij; 



d ^ - fi 
dr r 



= 0. 



(85) 
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It is necessary to solve this equation simultaneously with Eq. 
case takes the form 



which in this 



+ + = 0. 

77 



(86) 



5.3 General solution of evolution equation for nonsingular 
mode of perturbations 

Assuming 

^{r,r]) = R{r)e{r]) 
and dividing variables in Eq. 

kR" + ea^R = 0, 



we obtain ordinary differential equations 



ea Tj — D 



1 + K 



(1 + 3k)2 



e = o. 



(87) 
(88) 



To guarantee similarity of the solutions the sign of the division constant should 
be opposite to the barotrop coefficient 



e = -sgn(K). 
Solving Eq. dST]) we get 



n ^ • ^ a 

R — Ci sm —i=r + G2 cos —;=r. 



(89) 



(90) 



In order the function do not contain peculiarities at the origin of the coordinates 
\E'(r, 77) (condition (?)) it is necessary and sufficient in the solution (|90| C2 = 0, 
thus 



R(t) = C(a) sin 



where C(a)is an arbitrary constant. 
At K > Eq. ([HHl) has its solution 



(91) 



e(ry) = -^J,(a77) + -^n(«ry), 



(92) 



where Js{z) and Ys{z) are the Bessel functions of the first and second genuses 
respectively and 



1 

^=2 



5 + 3k 



1 + 3k 



1 

-, KG [-~l,+oo). 



(93) 
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At K < the Eq. solution is expressed by the way of the Bessel functions 
of imaginary argument Is{z) and Ks{z) (see, e.g. [10]): 

Q{V) - ^Isiav) + ^Ksiar^). (94) 

As far as the functions Js{z)/y/z and Ks(z)l\fz tend to infinity at z — ^ in 
order to get restricted ones at 77 —>■ it is necessary to set in the formulae ([92]) 
and (El 



= 0. (95) 

Thus nonsingular solution of evolutionary equation for perturbations can be 
written in the form 

LI (x 
^{r,T]) = — / C (a) sin— =rJs{ar])da, k > 0. (96) 
VV Jo \/\k\ 

At K < we obtain simultaneously 

LI (X 

^ir,i]) = — / C {a) sin— =rls{ar])da, k < 0. (97) 
VV Jo v|k| 
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